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ABSTRACT 



We present spectroscopic observations with Keck/LRIS of LBQS 0015+0239, 
a pair of quasars at z = 2.45 with a separation of A9 = 2.2 arcsec (projected 
linear distance of 17.8 hjQ kpc, for n„i = 0.3, A = 0.7). Lensing is an unlikely 
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interpretation for the images, since the spectra show significant differences in 
the N V and C IV emission fine profiles, and there is no luminous galaxy at the 
anticipated lens position. Rather, we interpret this pair as the highest rcdshift 
known example of a binary quasar. The redshift difi^crcncc of 661 ± 173 kms^^ 
between the two components is inconsistent with the lensing expectation, but is 
consistent with the line of sight velocity difference of a bound pair of galaxies. 
The spectra show associated, narrow metal absorption in highly ionized species 
of C and N, with a systemic velocity difference of 492 ± 6 kms"^. It is hkely that 
this absorption arises in a group or cluster of galaxies surrounding the quasar 
pair. After a thorough search, the Large Bright Quasar Survey (LBQS) is known 
to contain one gravitational lens and four probable binary pairs. The existence 
of close binaries is evidence that quasar activity can be triggered by galaxy in- 
teractions. We speculate that the close pair LBQS 0015+0239 is a rare example 
of the pre-merger state of two supermassive black holes. 

Subject headings: quasars: individual (LBQS 0015-1-0239) — quasars: emission 
lines — quasars: absorption lines — galaxies: interactions 

1. Introduction 

The rare situations of closely-separated quasars at the same rcdshift are of particular 
interest for observational cosmology. The paired images occur at a rate of ~ 1 in 500 quasars, 
and they are either the result of gravitational lensing or of physically distinct quasars in a 
binary system. Many of the pairs (and all of the triples and higher multiplicities) are the 
result of gravitational lensing by foreground mass distributions. Over 60 lens systems are 
known (see the review of Surdej & Claeskens 2002, and the constantly updated fist at the 
CASTLES web site http://cfa-www.harvard.edu/castlcs/). To be confirmed as a lens, the 
primary requirements are identical quasar redshifts, and a luminous galaxy at a plausible 
position and redshift to cause the observed image splitting and brightness ratio (Blandford 
& Narayan 1992; Schneider, Ehlers, & Falco 1992; Kochanck 1995). Almost all of these 
systems have separations smaller than A6 = 3 arcsec. At larger separations, there is a set of 
about 20 quasar pairs where the interpretation is not secure. A system is best confirmed as a 
binary if the spectral energy distributions are substantially different, including, for example, 
cases where one quasar in the pair is radio-loud and the other is radio-quiet. The binary 
population is important because it leads to the inference that quasar activity can be fuelled 
by galactic interactions and mergers (Djorgovski 1991; Peng et al. 1999; Kochanek, Falco, 
& Munoz 1999; Mortlock, Webster, & Francis 1999). 
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It is difficult to decide witfi any certainty wfietfier an individual pair is a lens or a 
binary. The presence of a luminous, and usually early type, galaxy along the line joining the 
two images is a good indication of lensing, since the chance probability of such alignment 
is very small. However, in some cases the available imaging is not deep enough to detect 
the putative lensing galaxy at all plausible redshifts, and the possibility of "dark" galaxies 
cannot be entirely ruled out (Hawkins 1997). Image separations of /S.6 = 1 to 3 arcsec are 
indicative of lensing but not prescriptive, since binaries in a final pre-merger state might 
have small projected separations. Conversely, the larger separations of = 3 to 10 arcsec 
need not be binaries, since galaxy lensing can get an assist from a surrounding cluster 
potential (Walsh, Carswell, & Weymann 1979; Lawrence et al. 1984; Wisotzki et al. 1993). 
Measurement of out-of-phasc but identical fractional flux variations between the two images 
would be sufficient evidence to declare a lens; but this type of data has not been acquired 
for any system where the lensing interpretation was not already secure. Another prediction 
of lensing is achromaticity. However, a clean test for spectral differences is complicated by 
the generally high degree of similarity of quasar spectra (Francis et al. 1992), the intrinsic 
variations in quasar continua and broad lines (Small, Sargent, & Steidel 1997), and the 
possibifity of differential extinction and reddening from a dusty lens (Falco et al. 1999). 
Spectral differences can also be caused by microlensing, where a stellar mass object can 
temporarily magnify the continuum of one image relative to the emission lines (Wisotzki et 
al. 1993; Lewis et al. 1998). 

The status of an individual pair can only be evaluated from the ensemble of its observed 
properties. Several recent papers have attempted to do this. Kochanek, Falco, & Munoz 
(1999) and Peng et al. (1999) both worked primarily with data gathered as part of the 
Cf A/ Arizona Space Telescope Lens Survey (CASTLES), and Mortlock, Webster, & Francis 
(1999) applied their own criteria to cataloged quasar pairs. These lists are in broad agree- 
ment, containing 13, 14, and 16 binaries respectively, after subtracting off the pairs where 
lenses have subsequently been detected. Two out of the 16 binaries listed by Mortlock et al. 
are noted as being possibly lensed. With Q 1208-1-1011, the luminosity and small separation 
are suggestive of lensing. Recent Hubble Space Telescope observations have failed to detect 
a lens, but the small separation and high contrast does not allow a sensitive limit (Lehar et 
al. 2000). Since there is no actual evidence that it is a binary, it is preferable to consider it 
as a probable lens. In the case of the wider separation pair Q 2345-1-007, lensing was intially 
indicated by the strikingly similar spectra (Steidel & Sargent 1991). However, recent Chan- 
dra observations find differences in the overall energy distributions that strongly indicate 
that this system is a binary (Green et al. 2002). The number of possible binaries increases 
from 15 to 19 if the objects discovered subsequent to the Mortlock et al. paper are included: 
CLASS B0827+525 (Koopmans et al. 2000), CTQ 839 (Morgan et al. 2000), the serendip- 
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itously discovered LBQS 0103-2753 (Junkkarinen et al. 2001), and LBQS 0015+0239 (this 
paper) . 

Another approach is to compare the statistical properties of the population of possible 
binaries to the properties of confirmed lens systems. The main caveat to this approach is 
the inhomogeneous nature of the sample of known pairs, which come from a wide variety of 
radio and optical surveys. The detection of one or both quasars as radio sources provides 
additional discriminatory power, since the radio-to-optical flux ratios of quasars vary over a 
wide range, and most quasars are radio-quiet (Hooper et al. 1996; Bischof & Becker 1997). 
Kochanek et al. (1999) have used the absence of pairs where both quasars are radio-loud, 
and the incidence of pairs where one quasar is radio-loud, to argue that essentially all the 
pairs where no lens is seen are in fact binaries. In addition, the distribution of redshifts, flux 
ratios and separations of the pairs is different from comparable distributions for conflrmed 
lenses (Mortlock et al. 1999). The spectral similarity of the quasar pairs is suggestive, but 
not compelling, evidence in the other direction (Peng et al. 1999; Mortlock et al. 1999). 

In this paper, we present observations of a small separation quasar pair that we interpret 
as a probable binary system. The new system is noteworthy because of the high redshift, 
the small angular separation, and the fact that it comes from a well-studied quasar catalog. 
Although the Large Bright Quasar Survey (LBQS) has been eclipsed in size by the 2dF and 
SDSS quasar surveys (Richards et al. 2001; Croom et al. 2001), it has been the subject of a 
large amount of follow-up activity, and its completeness has been calibrated with respect to 
the FIRST survey (Hewett, Foltz & Chaffee 2001). Most importantly, it has been thoroughly 
searched for small separation pairs (Hewett et al. 1998). In section 2, we present Keck/LRIS 
observations which resolve the components of this 2.2 arcsec separation pair. The spectra 
are noteworthy for the associated heavy element absorption that is seen in each component. 
In section 3, we present the evidence that LBQS 0015-1-0239 is a binary rather than a lens. 
Section 4 contains comments on the incidence of lenses and binaries in the LBQS, and section 
5 offers brief comments on the overall population of binary quasars. Conclusions follow in 
section 6. 



2. Observations 

All the LBQS quasars have been inspected visually on UK Schmidt Telescope plate 
material to provide a census of companion objects within a radius of lO'.'O (Hewett et al. 
1998). Companions with magnitudes brighter than Bj — 21.5 — 22, depending on the quality 
and depth of the Bj plates, have been identifled. However, in practice, the physical extent of 
the high density images of the relatively bright quasars on the plates precludes the detection 
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of faint companions closer than ~ 3'.'0 and LBQS 0015+0239 appeared as a single object 
on the original plate material. A supplementary program of broad band imaging at optical 
wavelengths, using a variety of telescopes, has produced CCD images of 230 LBQS quasars. 
The majority of the images were obtained under conditions of good seeing (< I'.'O), but are 
not particularly deep. The companion to LBQS 0015+0239 was identified in VRI images 
obtained in 1996 at the Im Jacobus Kapteyn Telescope and the image pair was flagged as a 
potential gravitational lens/binary quasar. 

The quasar pair LBQS 0015+0239 was observed on the night of 2001 July 21 with 
the 10m Keck 1 telescope and the Low Resolution Imaging Spectrograph (Rodgers & Mc- 
Carthy 1994). Observations were made in "Red Only' mode. The Red CCD is a 2048x2048 
SITe/Tektronix with 24/im pixels, giving 0.215 arcsec/pix averaged over the field of view. 
The 600 1/mm grating was centered just off its 5000A blaze at 5310A providing wavelength 
coverage of 4000-662lA with 1.28 A/pix dispersion and a resolution of 4.2A. The slit width 
was 1 arcsec and the slit was aligned with the two images. Three 15 minute exposures were 
obtained in approximately 0.7 arcsec seeing conditions, which was sufficient to resolve the 
two images. Flatfields were taken, a fiux standard was observed (BD+28°4211), and HgNeAr 
and CdZn arc lamps were obtained for both the science and calibration frames. 

The data was reduced using standard IRAF^ routines. The dual amp readout of the 
LRIS CCD necessitates the use of the Keck Observatory package wmkolris for overscan and 
bias subtraction. The fiatficlding required an extra step in the determination of the CCD 
response function because of the discontinuity introduced by the dual amp readout of the 
chip. Because the standard routines do not allow for discontinuities, we simply fit each half of 
the chip independently using the standard routines, then joined the two halves to obtain the 
complete response function. Before extracting the one-dimensional spectra, the three frames 
are shifted and median combined to eliminate cosmic rays. A variance weighted extraction 
is performed for each quasar component, including flux from ± 4 pixels to either side of the 
peak of the spatial proflle. To obtain a summed spatial proflle, the flnal combined image is 
corrected for geometrical distortion and pixels corresponding to AA 4600-6400 are summed 
along the dispersion axis using the IRAF task blkavg. All of the major absorption features 
seen in the combined spectrum are present in the individual 15 minute exposures. 

Spectra of the two components of LBQS 0015+0239 are shown in the panels of Figure 1. 
Prominent emission lines are labelled, and narrow absorption lines are marked along the 
continuum. The selection of absorption lines is based on software that uses a similar detection 



^The Image Reduction and Analysis Facilities package is distributed by NOAO, which is operated by 
AURA Inc. under contract to the National Science Foundation. 
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strategy to that used by the Quasar Absorption Line Key Project (Schneider et al. 1993). 

The basic features of the software are given in Petry, Impey, & Foltz (1998), and additional 
details can be found in Petry et al. (2002). The continuum is modelled with a spline fit 
that adequately follows the rapidly changing flux in the vicinity of the broad emission lines. 
Absorption features are selected using an interative technique that assumes unresolved lines 
and can deblend multiple Gaussian profiles in regions of complex absorption. 

Tables 1 and 2 list the significant absorption lines detected in the two components A 
and B of LBQS 0015+0239. Significance is defined in two ways: SLfu which refiects how 
well the absorber is fit by the software, and SLdet which refers to the detection limit of the 
data. SLfit is simply the ratio of the fitted equivalent width to its error, W/afu. SLdet is 
different in that its denominator, a^et-, is the la limiting equivalent width, which is a function 
of wavelength and is computed using the flux error. To obtain a reliable list of absorbers, we 
use a joint criterion for declaring significant lines: SLfu > 3 and SLdet > 4. Additionally, 
some lines in blended regions can have very low SL fu even though the detection limit of the 
data may be low, resulting in a higher SLdet- Therefore, we include any line with SLdet > 6. 
The RMS error in the wavelength calibration is 0.07A, although the absolute error is higher 
shortward of 4500A, where there are very few arc lines. Differential wavelength errors are 
smaller, but in practice the precision of measuring wavelength differences between the spectra 
is limited by the error on the fitted central wavelength of each feature, as can be seen in 
Table 1. 

The difference in signal-to-noise ratio between the two spectra refiects the fact that 
component A is nearly ten times brighter than component B. Each spectrum shows strong 
LycK, S IV, and C IV emission, with the red end of the Lya forest seen below 4200A. N V is 
strong in the wings of Lya in component A and virtually absent in component B. Strong and 
narrow associated absorption lines are seen eating into the Lya, N V, Si IV, and C IV emission 
features. The redshifts of the quasars were measured by cross-correlation with the SDSS 
composite, using the IRAF routine fxcor. The results are z — 2.4747 ±0.0043 for component 
A and z = 2.4831 ±0.0044 for component B. The difference of Av = 720 ±380 kms^^ differs 
from zero at the 2a level. Redshifts from cross-correlation are reliable because they use 
the full information in the spectra. However, the two components of LBQS 0015±0239 are 
more similar to each other than they are to the SDSS composite (as indicated by a higher 
correlation amplitude), so the most precise measure of the velocity difference comes from the 
cross-correlation between A and B. The cross-correlation result is Av — 661 ± 173 kms~^, 
which is formally different from zero at the Aa level. The amplitude of the cross-correlation 
is 0.92, indicating the high degree of spectral similarity of the two components. 

The major emission fines in each of the two components are stronger, by factors of 1.4 
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(Lya/N V), 2.3 (Si IV/0 IV), 2.0 (C IV), and 3.5 (He II) than the rest equivalent widths 
for the SDSS composite spectrum (Vanden Berk et al. 2001). The agreement between the 
two components is in all cases better than 35%. Peng et al. (1999) found that the spectral 
similarities of the binary pairs were significantly greater than would be expected for an 
equivalent sample of randomly selected quasars (but see also Mortlock et al., who reach 
a different conclusion for two of the LBQS pairs). The spectra of the two components of 
LBQS 0015+0239 raise again the question of how two quasars with an order of magnitude 
different luminosity could have such similar properties on parsec scales. 

The associated absorbers have a systematic velocity difference of 492 ±6 kms~-^ between 
the two components. This number is based on a velocity error weighted combination of the 
offsets between A and B of four different species, hsted in Table 3. The two components 
of C IV are combined into a single error-weighted offset before averaging with the other 
three species. Table 1 also hsts Lya absorbers seen at the short wavelength ends of the 
spectra, with 16 lines detected in component A, and 9 lines detected in component B. There 
are 5 coincident absorbers, where the wavelengths match within 2.5cr of the measurement 
errors. An additional 3 absorbers in B are resolved into 2 components each in A, where the 
average wavelength in A matches within the measurement errors to B. The non-detection of 
an absorber in A to pair to the remaining unmatched line in B is most probably due to the 
uncertain continuum fit on the wing of Lya emission in A. The strengths of these coincident 
Lya absorbers are consistent with being equal, given the different signal-to-noise ratios and 
equivalent width limits of the two spectra except for the two matches on the wing of Lya 
emission. This result is expected, since low column density hydrogen at this redshift has a 
coherence length much larger than the transverse separation of these sightlines (Dinshaw et 
al. 1994; Fang et al. 1996). 



3. Interpretation of LBQS 0015+0239 

As discussed in the introduction, it can be difficult to directly confirm an individual 
quasar pair as either a lens system or a binary. In this case, the systemic velocity difference 
of 661 ± 173 kms~^ between the two components differs from the lensing expectation at the 
3.8(T level. Additional evidence for a binary can come when one image of a quasar pair is 
a stronger radio source than the other (Kochanek et al. 1999). However, LBQS 0015+0239 
is not in the area of sky covered by FIRST, and NVSS only provides an upper limit to the 
combined flux density for both quasars of 2.5 mJy at 1.4 GHz (Condon ct al. 1998). The 
second way to search for achromaticity uses the optical spectra. The upper panel of Figure 2 
presents the spectra of the two components with an arbitrary scaling and offset to illustrate 
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the spectral similarity. The lower panel of Figure 2 shows the ratio of component A to 
component B where an interpolated continuum was used to replace the narrow associated 
absorption lines in each spectrum. The flux ratio is roughly a factor of ten with no evidence 
for any strong change with wavelength. We fit a straight line to the flux ratio, after truncating 
the data below 4400A, since the Lya forest is heavily absorbed and the flux ratio there will 
depend on resolution and the signal-to-noise ratio. The slope is zero to within the error 
of the fit. There are two systematic and significant departures from a flux ratio of ten. 
The first occurs just redward of Lya, refiecting the much stronger N V emission fine of 
component A. The other is centered on C IV emission, reflecting different line proflles for 
the two components. The N V difference between A and B supports the interpretation of 
the pair as a binary, since this difference is much larger than the temporal variations seen 
in individual quasars on timescales appropriate to a time delay under the lens hypothesis 
(SmaU et al. 1997). 

Another approach is to look for a lensing galaxy. No deep image has been taken of 
LBQS 0015+0239 under good enough seeing conditions to put a constraint on a lens. How- 
ever, the Keck spectra are well-resolved spatially, so they offer the opportunity of a mea- 
surement. The spectra were linearized on the final coadded image and then collapsed in 
wavelength to produce a ID spatial crosscut along a line joining components A and B. Fig- 
ure 3 shows the resulting spatial profile. A Gaussian fit to the two component centers in 
the crosscut yields an image separation of 10.2 pixels or 2.20 ± 0.03 arcsec. We choose a 
wavelength range of 4600-6400A for the crosscut, because it approximates the V passband 
and it avoids the edges of the CCD where the spectra fiare out shghtly due to defocussing. 
The bright component A has a magnitude of = 18.7 on the discovery plate material 
(Hcwctt, Foltz, & Chaffee 1995). Direct photometry yields V = 19.4 for component A and 
= 21.9 for B (Hewett & Foltz 2002, in preparation). Under the lensing hypothesis, the 
most probable lens redshift is z = 0.44 at a position between A and B but 90% of the 
way towards the faint component (C. Y. Peng 2002, private communication). We assume an 
early type galaxy with a singular isothermal sphere halo, since known lenses generally fit this 
description (Keeton, Kochanek, & Falco 1998). At z — 0.44, the velocity dispersion would 
be 247 kms~^. Scahng from a local Faber- Jackson relationship predicts a lens luminosity of 
1.6L*, which isV — 20.7 after (passively) de-evolving the fight to 2; = 0.44. 

A simple lens model therefore sets the expectation that the lens should be located 
between components A and B, and it should be much brighter than the faint component. 
Even without deep imaging, our spectroscopy is sufficient to rule this out. A Gaussian is 
actually a poor fit to the wings of the image profile in Figure 3. Instead, we omit the peak 
of the profile and fit a Gaussian to the lower slope and wing of component A in the direction 
away from component B, since a lens will not contribute significant light at that position. 
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Assuming symmetry, we fold this profile about the peak of A and use that shape for the 
wing of A in the direction of B. Component B is modelled by the same shape scaled by the 
ratio of peak fluxes. This procedure accounts for essentially all of the flux between the peaks 
of the two components, resulting in a strong limit on residual lens light. We can use the 
results to calculate the cross-contamination of the spectra presented in Figure 1: 6% of B 
is light from A and 0.03% of A is light from B. The upper limit on the amount of residual 
extended light at any position between A and B corresponds to V > 22.8. 

However, in a realistic lensing scenario, the lens would be 10 times closer to B than 
to A. The amount of residual extended light on or near the fainter component B gives a 
much more stringent upper bound oi V > 24.4. This constraint is only weakened in the 
unlikely scenario that the lens is compact enough to be unresolved and is positioned exactly 
in front of component B. However, the flux ratio in Figure 1 shows no achromatic signature 
from the admixture of red light from an early type galaxy, and in fact, the 4000A break 
would be seen at any lens redshift in the range 0.08 < z < 0.62. Therefore, any putative 
lens galaxy might be at least 30 times underluminous for its mass compared to galaxies in 
confirmed lens systems. Relative to the stellar mass-to-light ratio of an early type galaxy, 
this corresponds to a limits of {M/L)y > 200. However, since the probability distribution 
for the "most probable" lens redshift is fairly broad, these limits could be a factor of 3-5 
weaker. Since there is no other compelling evidence for massive dark galaxies, we prefer to 
interpret LBQS 0015-1-0239 as a binary quasar. 



4. Incidence of Binaries and Lenses in the LBQS 

The LBQS contains 1067 quasars, selected according to their spectral energy distribu- 
tions on objective prism plates from the UK Schmidt Telescope (Hewett, Foltz & Chaffee 
1995). The survey covers an apparent magnitude range of 16.0 < ttib < 18.85, a redshift 
range of 0.2 < z < 3.4, and is known to be no more than about 10% incomplete fainter than 
niB — 16.5 based on a comparison with the FIRST radio survey (Hewett, Foltz, & Chaffee 
2001). The LBQS has been subject to a specific search for companions within 10 arcsec of 
each quasar, down to level of ttib = 21.5 (Hewett et al. 1998). The ongoing 2dF and SDSS 
quasar surveys will eventually produce substantial new samples of close pairs. However, in 
each case, extensive follow-up observations are required. Both surveys use fibers which pre- 
clude the spectroscopic observation of pairs with separations under 20 arcsec. In addition, 
pairs less than 6 arcsec (for the 2dF) and 1.5 arcsec (for SDSS) have images that are merged 
on the original survey material. At present, the LBQS is the largest sample that has been 
closely inspected for pairs. 
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There are five quasar pairs with separations under 10 arcsec in the LBQS. One pair is a 
confirmed gravitational lens system, and the other four are probable binaries. Hewett et al. 
(1998) described the follow-up spectroscopy of plausible candidates for quasar companions. 
Quasar pairs with widely discrepant redshifts arc expected to occur by chance. Given the 
surface density of quasars down to rriB = 21.5, Hewett ct al. (1998) predicted 1.6 quasar 
pairs with unrelated redshifts in the LBQS, which is consistent with the one known (Hewett 
et al. 1994). LBQS 1009-0252 is confirmed as a lens system at z — 2.74, primarily due 
to the detection of a lens galaxy at an appropriate location between the two components 
(Claeskens et al. 2001). In addition, there is metal line absorption in the two spectra at 
z — 0.87, which is close to the lens redshift as predicted from a model and as estimated from 
the photometric properties. There is also a third image which is a quasar at z = 1.62, so this 
system qualifies as a second quasar pair with discrepant redshifts. The largest surveys for 
gravitational lenses arc radio surveys carried out at Jodrcll Bank and the VLA (Patnaik et 
al. 1992; Myers ct al. 1995; Phillips et al. 2001). With 18 lenses from 11,670 radio sources, 
roughly half of which are doubles, the presence of only one lens in the LBQS is consistent 
with expectations. 

LBQS 1429-0053 is a pair at z = 2.08 with an angular separation of A^^ = 5.1 arcsec 
(Hewett et al. 1989). Under the lensing hypothesis, this separation is most likely to be 
caused by a single massive galaxy with an assist from a surrounding cluster. However, 
Hewett et al. find no luminous object at the expected lens position, with a lower limit of 
2000 < (M/ L) < 4000, depending on cosmology. In addition, the spectra show a difference 
of ~50% in the relative strength of the Lya-N V complex, therefore LBQS 1429-008 is a 
probably binary. LBQS 2153-2056 is a pair at z = 1.85 with an angular separation of 
A^^ = 7.8 arcsec (Hewett ct al. 1998). Available imaging is not deep enough to place 
a sensitive limit on a deflector, but the fact that only one member of the pair is a BAL 
quasar argues against lensing unless the BAL clouds arc extremely small (see Mortlock et 
al. 1999). LBQS 0103-2753 is by far the smallest separation binary known (Junkkarinen et 
al. 2001), Bit z — 0.85 with a separation of A^ = 0.3 arcsec. The binary interpretation is 
secure because the spectra are very different (and one quasar has BAL features). Under the 
lensing hypothesis, the time delay is ~10 days, and spectral variations sufficient to explain 
the differences between two components' spectra are never observed in quasars on such a 
short timescale. As we have argued in this paper, LBQS 0015-1-0239 is a strong candidate 
for a fourth binary quasar in the LBQS. 

The statistical properties of the LBQS pairs provide additional evidence that they are 
in fact binaries. Three out of four have flux ratios of ten or greater, versus only 3 out 
of 24 two- image lenses with data tabulated by CASTLES. Lack of spectral similarity pro- 
vides an additional argument. The primary of the serendipitously discovered close pair 
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LBQS 0103-2753 was selected as a BAL quasar, but the companion has a very different, and 
unexceptional, spectrum with no strong absorption. Mortlock et al. (1999) used principal 
compnents analysis to show that LBQS 1429-0053 and LBQS 2153-2056 are as dissimilar 
as any quasar pairs randomly selectly from the LBQS. This would not be expected if either 
pair were lensed. 

The pair that is the subject of this paper is unusual in that both quasars show associ- 
ated absorption. This is noteworthy, given that fewer than 10% of LBQS quasars have this 
attribute. Associated absorbers are usually defined as metal systems within Af = ±5000 
kms~^ of the emission redshift. Such absorption might be intrinsic and related to the broad 
absorption line (BAL) phenomenon, or extrinsic and caused by galaxies in a cluster con- 
taining the quasar (Foltz et al. 1988). Evidence for the latter interpretation comes from 
abundance studies. Nine out of ten associated systems have {N/ C) ratios below solar abun- 
dance, while BAL regions have (N/C) ratios of 9-10 times solar abundance (Turnshek et al. 
1996; Franceschini & Gratton 1997). This new pair also supports the extrinsic explanation 
of associated absorption, as the redshift difference of 492 ±6 kms~^ is consistent with a group 
or cluster environment, and the fact that the pair is a binary system shows that there are at 
least two massive galaxies in this region of space. 

To date, four binary quasars and a single gravitational lens have been discovered in the 
LBQS. The incidence of binary pairs is difficult to predict a priori, but it is worth asking 
what is expected for the rate of lensing in a survey like the LBQS. Maoz et al. (1993) 
conducted the largest optical lens survey that was sensitive to the angular separations of 
essentially all known lens systems. The survey found 5 lenses out of 502 quasars surveyed 
(Maoz gave a range of 3-6, so this count includes 0957-1-561 but excludes the unconfirmed 
candidate UM 1120-1-0154), for a lensing rate of 1%. For the purposes of this paper, we will 
not carry out a full analysis of the predicted lensing rate of the LBQS, but we can easily 
give an estimate relative to the Maoz HST snapshot survey. 

The predicted number of multiply imaged quasars for any sample is found by calculating 
the probability of each quasar being lensed and summing over these probabilities for the 
entire sample. The probability that an individual quasar is lensed is the product of its 
optical depth to lensing and its magnification bias (Turner, Ostriker, & Gott 1984; Fukugita 
& Turner 1991). Pei (1995) has published maps of magnification bias in redshift and apparent 
blue magnitude, assuming a double power law model for the quasar luminosity function and 
a lens population that is composed of compact objects, individual galaxies, and clusters. 
The shape of the magnification bias contours is insensitive to quite varied models of the 
lensing population. Taking the 50% and 90% bounds in redshift and magnitude for the two 
samples, it turns out that the mean magnification bias is low for both the LBQS and the HST 
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snapshot survey, just under 1%. The most substantial difference between the two samples 
is the mean redshift, about 1.2 for the LBQS and 2.0 for the HST snapshot survey. For a 
A-dominated cosmology, the lensing optical depth scales as t ^ (Bahcall et al. 1992). 
Thus, the lensing rate should be 4-5 times lower for the LBQS than for the HST snapshot 
survey, consistent with the single lens that has so far been identified. 

Despite this argument, it is admitted that the predictions of lensing rates are uncer- 
tain, and the purely observational limits on additional lenses in the LBQS are poor. The 
companion search of Hewett et al. (1998) was only sensitive to separations greater than 3 
arcscc (although equal brightness doubles as close as 2 arcscc could have been identified). 
The LBQS search was therefore insensitive to about 90% of the lens image separations from 
a radio survey such as CLASS. The presence of a few more small separation lenses in the 
LBQS cannot be ruled out. 

5. The Population of Binary Quasars 

The LBQS contains four binary pairs with projected separations less than lOOhjQ kpc 
(plus the previously mentioned lens which is, of course, a single quasar). Since the newly dis- 
covered LBQS 0015-1-0239 has smaller angular separation than the 3-10 arcsec completeness 
range of the search by Hewett et al. (1998), and since the 0.3 arcsec pair LBQS 0103-2753 
was discovered serendipitously (Junkkarinen et al. 2001), it is worth asking how many more 
small separation binaries might be lurking in the LBQS. Maoz ct al. (1993) performed a 
snapshot survey of 498 quasars with z > 1 and found no binaries with separations below 
3 arcsec. They could have detected binaries as close as 0.3 arcsec, and with fiux ratios of 
at least 4:1. Thus, the single close LBQS binary is consistent with the null result of Maoz 
et al., and Poisson statistics suggest that there are unlikely to be more than two additional 
very close binaries in the LBQS. 

Over the order of magnitude larger range of scales from 10 to 100 arcsec, there are no 
LBQS pairs with similar rcdshifts. Unfortunately, this null result is not a strong constraint, 
since the deep companion search of Hewett et al. (1998) applied over this angular range 
would have been been sensitive to 30 times as many quasars. As noted originally by Djor- 
govski (1991), the high incidence of close pairs is strong circumstantial evidence that quasar 
activity can be triggered by direct interactions. Extrapolating the observed quasar corre- 
lation function down to these small physical scales leads to a prediction of ~ 0.01 quasar 
pairs within lOOh^Q kpc, compared to the four observed. A similar scaling of the correlation 
function predicts 20 times as many pairs between 10 and 100 arcsec as between 3 and 10 
arcsec, or ~40 wide pairs. 
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We adopt the pair fraction of the LBQS (4/1055) as the probabihty that a quasar is a 
member of a binary, Punary = 4 x 10~^. The LBQS probabihty only apphes out to 10 arcsec, 
or a transverse separation of ~ lOO/ifo^ kpc. If the interaction-induced quasar activity occurs 
in the advanced stages of a merger process, then we anticpate that the probabihty summed 
to much larger separation scales will not be much larger (Barnes & Hernquist 1996). If 
the radio-loud fraction of optically selected quasars is frad and we assume that binarism 
does not affect radio properties, then the relative numbers of O^, O^i?, and O^i?^ binaries 
(ie. those where both, one, and neither quasar are radio-loud) are in the ratio (1 — fradf '■ 
'^fradi^ — frad) '■ /^arf (Kochanek et al. 1999). Updating the Mortlock et al. list (see comments 
in the Introduction) gives 12 pairs. Taking f^ad — 0.1, as is appropriate for FIRST or 
NVSS detections above 1-3 mJy, the predictions arc 0.1 O^B? pair, which is consistent with 
none seen, and 2.2 O^i?, which is consistent with two seen. 

Recent radio lens surveys are an order of magnitude larger than the LBQS. Also, the 
surface density of radio sources is low enough that it is much more feasible to search for 
pairs out to larger angular separations. So far, only one radio-loud binary {O^R^) has been 
found among 15,000 CLASS sources (Koopmans et al. 2000), and none have been found 
among 9100 FIRST sources (Ofek et al. 2002). The only other known O^R^ binary is 
MG 0023-1-171 (Hewitt et al. 1987). It is entirely plausible that several more O^R^ binaries 
remain to be found since the radio morphologies of each component may be quite different, 
thus eluding the selection designed to catch lens candidates. Using the LBQS estimate of 
Pbinary: the total number of O'^R binaries present in the CLASS and FIRST surveys out to 
= 10 arcsec should be ~57 and ~35, respectively. The number where both components 
are radio-loud is a factor of frad smaller, or ~3 and ~4, respectively. The numbers should 
also be increased to account for the larger angular scales of the CLASS (15 arcsec) and 
FIRST (30 arcsec) surveys, assuming that the association probability scales according to 
the quasar-galaxy correlation function. The smaller size of the FIRST survey is outweighed 
by the deeper radio limit and the larger angular search radius. For the FIRST survey, we 
predict ~15 O'^R binaries. We predict ~100 and ~400 binaries in the 2dF and SDSS optical 
quasar surveys, but note that substantial spectroscopic foUowup will be needed to identify 
the second (fainter) components in each survey. 

LBQS 0015-1-0239 is the third smallest separation binary quasar known. Kochanek et al. 
(1999) have shown that number of binaries is consistent with the quasar-galaxy correlation 
function at small scales. Mortlock et al. (1999) deprojected the observed separations of 
binary pairs to yield a broad distribution peaking at ~ 50h^Q kpc, a declining tail out to large 
separations, and a possible deficit at small separations. Despite the lack of comprehensive 
follow-up of wide pairs from radio surveys, the decline at large separations appears to be 
real and is indicative of an interaction trigger for the nuclear activity. LBQS 0015-1-0239 and 
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the other two binary quasars discovered since Mortlock et al. neatly fill in their "hole," so 
the separation distribution may be flat below bOh^Q kpc. Given this type of distribution, 
the 0.3 arcsec LBQS pair discovered by Junkkarinen et al. (2001) is far more likely to be 
a real 2-3 hj^ binary than a chance projection of a ~ AOhjQ kpc binary. LBQS quasars 
are typically found in massive galaxies (Hooper, Impey & Foltz 1997), where the dynamical 
decay timescale of the binary orbits is close to a Hubble time and much longer than active 
timescale of quasars accreting at near the Eddington hmit (Binney & Tremaine 1987). Thus, 
LBQS 0015+0239 may have been triggered by the much faster settling process with the 
merged halo (Barnes & Hernquist 1996). Deep imaging of this pair, and many of the other 
binaries, should reveal heavily distorted host morphologies. 



6. Conclusions 

We have carried out Keck/LRIS spectroscopy of the 2.2 arcsec separation quasar pair 
LBQS 0015+0239. These data, plus a consideration of the other pairs in the LBQS, leads 
to the following conclusions: 

1. There are three pieces of evidence that support the interpretation that the pair 
is a binary system. The velocity difference between the components is inconsistent with 
the lensing expectation but consistent with a bound pair of massive galaxies. There are 
substantial differences in N V, and to a lesser extent C IV, emission line strengths between 
the two components. A spatial cut along the line connecting the pair in the equivalent of 
the V band leads to the conclusion that any lensing galaxy must be at least thirty times less 
luminous than an elliptical galaxy that could account for the image separation. 

2. Narrow associated absorption due to metal lines is seen in each component, with a 
velocity difference of 492 ± 6 kms~^. Since the associated absorption phenomenon is only 
seen in about 10% of LBQS spectra, its presence in these two quasars supports the idea that 
the absorption arises in a group or cluster environment that includes the binary pair. 

3. There are four probable binary pairs in the LBQS, so the incidence rate is 4 x 10~^, 
several times higher than the lensing rate. The large radio surveys being used to select 
gravitational lenses probably contain 15-20 unrecognized binary pairs, where only one of the 
members is radio-loud. LBQS 0015+0239 has a projected transverse separation of only 17.8 

kpc and is the highest redshift binary quasar known. It may represent the short-lived 
phase prior to the coalescence of two supermassive black holes, in which case deep imaging 
of the pair should reveal a disturbed morphology for the host galaxies. 
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Fig. 1. — Keck/LRIS spectra of the quasar pair LBQS 0015+0239, with 4.2A resolution. 
Cross correlation with SDSS quasar composite gives redshifts of z = 2.4747 for component 
A and z = 2.4831 for component B. The dotted line is the fitted continuum, prominent 
emission features are labelled, and all significant absorption lines are marked with every 
other one numbered. The lower curve shows the 3a flux error. The flux is in units of 10"^^ 
ergs s~^ cm~^ A~^. 
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Fig. 2. — The upper panel shows the spectra of LBQS 0015+0239 components A and B 
scaled and offset by an arbitrary amount. The lower panel shows the ratio of component 
A to component B after interpolation over the narrow associated absorption features. The 
dotted line is plotted at the approximate mean flux ratio. 




Fig. 3. — The upper panel shows the spatial cross-cut of LBQS 0015+0239 over the spec- 
tral range 4600-6400A. The fitted profile (dashed line) is scaled from the left-hand side of 
component A, where any potential light from a lensing galaxy will be negligible. The lower 
panel shows the residual counts after subtraction of the two fitted components. In both 
panels the counts are in units of 10^, and the vertical dotted hues show the centers of the 
two components. 
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Table 1. LBQS 0015+0239 A 
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Table 1 — Continued 
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5, 


.02 


43 


5889. 


,74 


± 


0. 


,36 


0, 


,520 


± 


0, 


,070 


4.20 ± 0.00 


1, 


.03 


7, 


.00 


8, 


.68 


44 


5895. 


,53 


± 


0. 


,60 


0. 


,290 


± 


0. 


,070 


4.20 ± 0.00 


1, 


.03 


4, 


.30 


4, 


.85 



'^Wavelengths are vacuum heliocentric. 

*^Lincs with a value of 4.20 ± 0.00 were fit with the FWHM set to the minimum 
allowed value. 

*^ Significance of the line defined as EW/afu, where afu is the error in the measured 
equivalent width. 

"^Significance of the line defined as EW/adet, where aa^t is the la limiting equiva- 
lent width. 

*^This line has been identified with an associated absorption system and is included 
in Table 3. 
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Table 2. LBQS 0015+0239 B 



Line 


Central A'' 


Equivalent Width 






(A) 


(A) 


(A) 



1 


4055, 


OH' 

.07 


1 


0, 


.79 


2. 


.730 


1 


0, 


o oo 

.820 


5, 


.56 


1 


o 

2, 


OO 

.02 


1 


o^ 

.07 


3. 


OO 

.30 


5. 


1 1 

.11 


2 


4062, 


.30 


1 


0, 


.65 


2. 


,390 


1 





n r\ 

.560 


4, 


oo 

.20 


1 


0, 


oo 

.00 


1, 


o^ 

.07 


A 

4. 


oo 

,20 


A 

4. 


A O 

,42 


3 


4097, 


oo 

.22 


1 


1, 


O 1 

.21 


2. 


,890 


1 


0, 


.810 


8, 


.87 


1 


2, 


oo 

.92 


i 

1, 


.16 


3, 


,60 


6. 


/I o 

,40 


4 


A "1 OO 

4122. 


.33 


1 

± 


0, 


.97 


1. 


,590 


1 


0, 


.530 


6, 


.01 


1 


2, 


.37 


0, 


.35 


3. 


oo 

,00 


4. 


A n 

,46 


5 


4133, 


A O 

.43 


± 


0, 


.37 


4. 


,310 


± 


0, 


.500 


6, 


.68 


± 


0, 


.94 


1, 


.10 


8. 


,60 


13, 


.52 


6 


4147, 


.43 


± 


0, 


.97 


2. 


,700 


± 


0, 


.580 


9, 


.62 


± 


2, 


.55 


0, 


.92 


4. 


,70 


9. 


,13 


7 


4191, 


.85 


1 


0, 


OO 

.32 


o 

2. 


.740 


1 


o 

0, 


.310 


5, 


f7 A 

.74 


1 


o 

0, 


.77 


1 


A O 

.49 


8. 


OO 

.80 


1 o 

12, 


r o 

.59 


8 


4200, 


.63 


1 


0, 


.55 


1. 


.170 


1 

± 


0, 


O 1 o 

.210 


4, 


.20 


1 

± 


0, 


AA 

.00 


1 


y1 A 

.49 


5. 


cr A 

.50 


D. 


OT 

.37 


9 


4204, 


.79 


± 


0, 


.50 


1. 


.130 


± 


0, 


.210 


4, 


.20 


± 


0, 


.00 


1 


.49 


5. 


.40 


6. 


.92 


10^ 


4221, 


.22 


± 


0, 


.58 


2. 


.910 


± 


0, 


.350 


4, 


.20 


± 


0, 


.00 


0, 


.88 


8. 


.40 


24, 


.42 


11 


4224, 


.33 


± 


0, 


.41 


1. 


,930 


± 


0, 


.440 


4, 


.20 


± 


0, 


.00 


0, 


.88 


4. 


,40 


16, 


.68 


12<= 


4301, 


.56 


± 


0, 


.26 


2. 


,150 


± 


0, 


.260 


4, 


.34 


± 


0, 


.61 


0, 


.64 


8. 


,20 


9. 


,91 


13*^ 


4315, 


.10 


± 


0, 


.51 


1. 


,180 


± 


0, 


.300 


4, 


.20 


± 


1, 


.29 


1, 


.56 


3, 


,90 


4, 


,99 


14e 


4839, 


.88 


± 


0, 


.62 


1. 


,540 


± 


0, 


.370 


5, 


.44 


± 


1, 


.56 


1, 


.21 


4. 


,20 


6, 


,03 


15^ 


4870, 


.19 


± 


0, 


.35 


1. 


,560 


± 


0, 


.210 


4, 


.20 


± 


0, 


.00 


0, 


.77 


7. 


,40 


7. 


,17 


16^ 


5375, 


.40 


± 


0, 


.13 


6. 


.370 


± 


0, 


.250 


7, 


.30 


± 


0, 


.35 


2, 


.42 


25, 


.50 


50, 


.60 


17^ 


5384, 


.23 


± 


0, 


.13 


4. 


.080 


± 


0, 


.220 


5, 


.34 


± 


0, 


.31 


2, 


.42 


18, 


.60 


33 


.04 


18 


5712, 


.25 


± 


1 


.18 


1. 


.170 


± 


0. 


.390 


7, 


.29 


± 


2. 


.83 


0, 


.40 


3. 


.00 


4. 


.86 


19 


6330 


.57 


± 


0, 


.59 


1. 


,600 


± 


0. 


,380 


4, 


.20 


± 


0. 


,00 


0, 


.44 


4. 


,20 


4. 


,55 


20 


6343, 


.12 


± 


0, 


.48 


1. 


,670 


± 


0. 


,320 


4, 


.20 


± 


0. 


,00 


1, 


.77 


5. 


,30 


5. 


,16 


21 


6585, 


.49 


± 


0, 


.88 


1. 


,120 


± 


0. 


.380 


5, 


.39 


± 


2. 


,10 


0, 


.09 


3. 


,00 


4. 


,01 



^Wavelengths are vacuum heliocentric. 

^Lines with a value of 4.20 ± 0.00 were fit with the FWHM set to the minimum 
allowed value. 

'^Significance of the line defined as EW/afu, where a fa is the error in the mea- 
sured equivalent width. 

'^Significance of the line defined as EW/adet, where cr^et is the la limiting equiv- 
alent width. 

^This line has been identified with an associated absorption system and is in- 
cluded in Table 3. 
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Table 3. Associated Absorption 





Component A 


Component B 






Line 


Central A 


Redshift 


Central A 


Redshift 


Identification 


A V 




(A) 




(A) 






(km s~^) 


1 


4228.10 ± 0.02 


2.4780 


4221.22 ± 0.58 


2.4723 


H I 1216 


488.2 ± 41.2 


2 


4306.59 ± 0.12 


2.4764 


4301.56 ± 0.26 


2.4723 


N V 1238 


350.3 ± 20.0 


3 


a 




4315.10 ± 0.51 


2.4721 


N V 1242 




4 


4846.76 ± 0.67 


2.4775 


4839.88 ± 0.62 


2.4726 


Si IV 1393 


425.9 ± 56.5 


5 


b 




4870.19 ± 0.35 


2.4718 


Si IV 1402 




6 


5384.81 ± 0.06 


2.4811 


5375.40 ± 0.13 


2.4720 


C IV 1548 


524.3 ± 8.0 


7 


5392.90 ± 0.08 


2.4776 


5384.23 ± 0.13 


2.4720 


C IV 1550 


482.3 ± 8.5 



^ Both components of the N V doublet are detected in absorption in Component B. However, 
because N V emission is much stronger in Component A, ambiguities in continuum fitting have 
made it difficult to detect N V 1242 in absorption for this object. 

^ This absorption feature was detected in the spectrum of the A component, but it did not 
meet the significance criteria and is not included in Table 1. 



